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ABSTRACT 

Protons gain energy in short range collisions with heavier dark matter particles (DMPs) of comparable 
velocity dispersion. We examine the conditions under which the heating of baryons by scattering off DMPs can 
offset radiative cooling in the cores of galaxy clusters. Collisions with a constant cross section, a^^ , independent 
of the relative velocity of the colliding particles, cannot produce stable thermal balance. In this case, avoiding 
an unrealistic increase of the central temperatures yields the upper bound a^^ < 10^^^ cm^{mx/mp), where 
TTix and rrip are the DMP and proton mass, respectively. A stable balance, however, can be achieved for a 
power law dependence on the relative velocity, V , of the form, a^^ oc V"" with a < —3. An advantage of this 
heating mechanism is that it preserves the metal gradients observed in clusters. 

Subject headings: cosmology: dark matter - galaxies - galaxy clusters: cooling flows 



1. Introduction 

While the existence of "dark" matter particles (DMPs) 
is one of the main pillars of the standard cosmological 
paradigm, the nature of these particles and the way 
they interact with the ordinary matter are still un- 
known. Presently the weakly interacting massive particle 
(WIMP) is the most popular candidate for the role of the 
DMP, but as long as it remains undetected other can- 
didates present viable alternatives. Strongly interacting 
heavy particles are produced non-thermally in some cos- 
mological scenarios (e.g. Griest & Kamionkowski 1990; 
Chung, Kolb & Riotto 1998). Underground experiments, 
which are primarily designed to search for WIMPs, are 
insensitive to particles whose interaction cross-section 
with baryons, a^^, is above ^ 10^''*' cm^. Since the 
flux of such particles is heavily diluted by the ground or 
even by the atmosphere (Albuquerque & Baudis 2003), 
the most rigorous constraints are based on satellite ex- 
periments (Wandelt et al. 2000) and observations of 
astrophysical systems (Starkman et al. 1990; Chen et al. 
2002) Current constraints allow for super-heavy particles 
with masses above 10^ Gev with a wide range of cross 
sections with ordinary matter. On the low side, DMPs 



with mx < 0.4 Gev are still admissible. 

In this paper we examine the role of short range inter- 
actions between the DMPs and protons in galaxy clus- 
ters. As particles fall into the halos, the gravitational 
virial theorem implies that baryons and DMPs acquire 
similar velocities, V . Since the kinetic energy of a parti- 
cle of mass m is mT^^/2, the direction of coUisional energy 
transfer between baryons and DMPs depends mostly on 
their masses. For > rrip and large enough coUisional 
cross section, heating by DM overcomes radiative cooling 
and increases the temperature. If this were the case for 
galaxy clusters, then we would expect to find the largest 
temperatures at their centers, where the density and col- 
lision rates are highest. However, in most clusters the 
observed temperature decreases towards the center. This 
allows us to place a new upper limit on a^^ . For veloc- 
ity independent scattering the new limit is two orders of 
magnitude lower than the previous limit obtained from 
stability analysis of galactic disks and DM halos (7). 

Qin & Wu (2001) have argued that for :§> rup 
and CTxp ~ 10~^^(mx/wp) cm^, which is similar to the 
range of strong interaction, energy transfer from DM to 
baryons can explain the deficit of cold gas (T < 10^ K) in 
the cores of clusters of galaxies, as suggested by Chandra 



and XMM-Ncwton observations (Kaastra 2001, Peterson 
2001, Tamura 2001). Here we elaborate on the work of 
Qin and Wu (2001) in several ways. We use a more ac- 
curate description of the radiative cooUng and heating 
rates and consider the consequences of a velocity depen- 
dent cross section. Further, we account for the increase 
in the velocity dispersion of the dark matter as it con- 
tracts while transferring energy into the protons. We 
show that unless the cross-section strongly varies with 
energy, a state of thermal equilibrium between heating 
by the DMPs and radiative cooling is unstable. 

The paper is organized as follows. In §2 we write down 
the relevant equations. In §3 we check the requirements 
for the stable thermal balance between radiative cooling 
and heating by DMPs. In §4 we derive analytic upper 
limits on the cross sections. We present numerical solu- 
tions in §5 and conclude with a discussion in §6. 

2. The equations 

We begin by writing down the equations governing the 

temperature and density of a gas element subject to ra- 
diative cooling and heating by collisions with dark matter 
particles. As the self-gravitating dark matter particles 
transfer energy to the baryons, they get closer to the 
center and become hotter. Although the mass dark mat- 
ter particle is much larger than that of the proton, the 
energy density is not and therefore we must account for 
the increase in the velocity dispersion of the dark matter 
particles. We denote the gas heating and cooling rates, 
both per unit volume, by _ffx and L, respectively. The 
corresponding rates per unit mass are H^/rit and L/nt 
where rit ~ 2.2np is the total number density of all parti- 
cles in a fully ionized primordial plasma. The energy bal- 
ance equation is most conveniently expressed in terms of 
the entropy. We define an "entropy" as 5 = k^T/rie^/^ 
where rio « 1.1 6np is the electron number density in 
cm~'^ and k^T is in kev. The entropy obeys. 
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The advantage of this form of the energy balance equa- 
tion is that the adiabatic heating/cooling is readily in- 
cluded. 

As will be shown below, in the first approximation 
the heating rate i?x depends on the dark matter prop- 
erties only through the energy density (per unit volume) 
E-^ = {l/2)n^m^V^ where is the number density of 
dark matter particles, and is their velocity dispersion. 
According to the virial theorem = — i^totai- There- 
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There are three time-scales in tlie probk;m, t^, 



(2) 

and 



i^jj which correspond, respectively, to the radiative cool- 
ing, heating of gas by collisions with DM particles, t^^, 
and kinetic energy gain of dark matter. These time scale 
are given by 



t = 



2L 

E^ E, 
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The virial theorem implies that after the collapse the 
energy densities of DM and gas relate as their mass den- 
sities, -EgasZ-Ex ~ ~ 0.15. However, the energy 
acquired by protons during the collapse is very quickly 
shared with electrons and helium ions so that the velocity 
dispersion of these species is fixed by energy equiparti- 
tion. Therefore, initially we have nprn-pV^ / n^myiV^ = 
{np/nt)fg^, ~ 0.07. Subsequent effects of DM-baryon in- 
teraction are then determined by the equations (1) and 
(2). The cross section must satisfy tlie condition 



(6) 



otherwise the temperature would increase at the cen- 
tral regions, contrary to observations. The last condition 
is sustained by hydrodynamical simulations which show 
that as clusters form their temperature profiles increase 
towards the center (e.g. Navarro, Frenk & White 1995). 
Therefore, if the heating by DMPs is stronger than cool- 
ing then the temperature will always be highest at the 
center. 

If the mass and typical kinetic energy of a DMP are 
much larger than those of a proton, rrix ^ m-p and 
fnyV^ ^ nipVp, then the DMP recoil in a single scat- 
tering with proton can be neglected. In this case, the gas 
heating rate per unit volume is (Qin & Wu 2001) 

27ry3y3 

X j d cos e {s"^ +v'^ -u'^)s, (7) 

where we have assumed Maxwellian energy distribu- 
tion for the DMPs and protons, and Vp are, respec- 
tively, the velocity dispersions of DMPs and protons, and 
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COS 9 = + — s'^)/2uv. Assuming a power law de- 
pendence , cTxp = o"o('S/Vb)", of the cross-section on the 
relative speed s = |Vi — Vp| transforms eq. 7 into 

2nV^V^ Jo 

X dcOS0(s2+u2_y2)^l+a_ ^g^ 

For integer values of a the above integrals can be solved 
analytically. However, the solutions are quite cumber- 
some, so it is more convenient to use an approximation 

« 6 • 2«mpnpn,aoFo""K'+"(l + ^) > (9) 

In the central regions of clusters Vx < Vp, which allows 
further simplification of eq. 9 

ifx » 6 • 2«mpnpnxCTo'5^o""K^^p, (10) 

We write the cooling rate as (Tozzi & Norman 2001) 

as 

L = nine [Ci(fcBT)« + C^ik^Tf + C3] (11) 

in units of 10~^^erg cm~^ s~^. The number density of 
ions ni ~ l.OSnp and electrons ng « 1.16np are in cm^"^. 
For metallicity of one third solar, a = —1.7, /3 = 0.5, 
Ci = 8.6 X 10-3, C2 = 5.8 X 10-2, and C3 = 6.4 x 10-^ 
and k^T is in kev. 

3. Thermal stability analysis 

Because the term dE^/dt never vanishes (see equation 
2), a strict thermal equilibrium cannot hold. However, 
the time scale, t^^, for significant changes in is typi- 
cally a factor of 1/ f^^^ « 6 longer than the time scale, t^j, , 
for heating the baryons. Therefore, it is meaningful to 
address the question of whether a quasi-equilibrium with 
AS/dt K, (see equation 1) is thermally stable over time- 
scales shorter than t^^ . We express H-x_ and L in terms of 
T and the pressure p = ntfeeT. We take Hy, = ApT"-/"^'^ 
and L = BT^~^p^ where A and B represents all quanti- 
ties that are independent of p and T. Further, we write 
InS = InT^/^/p + const. Equation (1) is then 

3^dln(TV3/p) ^ ^^^„/2-l/2 _ 5p2y/3-2 _ (12) 

Let the r.h.s vanish for T = Tq at p = po and consider 
variations in T at constant pressure p = po- Substi- 
tuting T = To + ST, the last equation gives (5/2)pT = 



(a/2+3/2-/3)ApT''/2-i/2 ^^.^^ ^^^.^^.j. gj. Therefore, 
the system is stable for a < 2/3 - 3. For T ft! lO^K and 
metallicity of one third solar, « 0.5, so that a < —2. 
For clusters with lower temperatures and higher metal- 
licity P « —0.3, and we get a < —3. Therefore, the 
dependence of the cross section must be steep in order to 
obtain thermal stability over time scales < t^^ . 

4. Upper limits on cross sections 

To allow the temperature at the core to fall the cooling 
must overcome the heating, L > H. Since both L and 
H are not constant, the value of a^^ which gives the 
thermal balance changes over time. However, unless the 
condition L > H \s not satisfied from the start, it will 
never be satisfied later. Thus we may derive the upper 
limit on the cross-section by using the initial conditions 
of clusters. 

For a > —2.3 the stringiest limit is obtained from 
the hottest clusters where the cooling is dominated by 
bremsstrahlung and the cooling function is well approx- 
imated by 

i = 2.9 x 10-2^n^T°-^ erg s-^ cm^ , (13) 
where T is in K. The limit on a is obtained from the 

xp 

condition (6) using the expressions (10) and (13). Taking 
UpTripVp / UxTrixV^ ~ 0.07, the condition (6) yields^ 

a < 2.7X 10-2 W .2-° V7 Ji-) " 

\mpj VlO^Ky \lm/sj 

(14) 

For energy independent scattering (a = 0) we get 

.,.<2.7xI0-w(jX_)"(;=^), (15) 

which for the hottest clusters, such as RX J1347.5-1145 
with T > 1.7 • 10^ K (Allen, Scmidt & Fabian 2002), 
yields 

a^^ < 1.5 X 10-25cm2 (^^'j . (16) 

For a < -2.3 it is the coldest clusters (T - 3 • 10^ K) 
which give the stringiest limit. There the eq. 13 is no 
longer a very good approximation and the full cooling 
function has to be used. Therefore we present the limit 
graphically in fig. 2. 



The upper limits on the cross sections do not depend on the actual 
value of Up since for a given f^^^ both the initial cooling and heating 
rate are proportional to rip. 
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The above calculations ignore several additional heat- 
ing processes such as AGN, thermal conduction and SN 
feedback. However, if these are important they will make 
our limits even more restrictive. 

5. Numerical solutions 

We show now results of numerical integration of the 
equations (1) and (2), for various values of initial tem- 
peratures T using the full expression for the cooling rate 
given by (11). Our goal is to confirm the upper limit de- 
rived in the previous section and to demonstrate that a 
reasonable temperature evolution is obtained for energy 
dependent cross sections. 

We would like to integrate the equations with initial 

conditions describing the physical state in the cool cores 
of clusters before cooling and heating become important. 
The initial temperature in the core can be extrapolated 
from the outer regions assuming an initial isothermal 
profile. We also need the initial number density, rip, 
in the core. This is more tricky to obtain as the ob- 
served rip is likely to be substantially different from its 
initial value and from the measured value in the outer 
regions. To illustrate the effect of varying rip we exper- 
iment here with two choices for rip. We use three val- 
ues for the initial temperatures, T = 16, 8 and 3 kev. 
Candidates for clusters with these temperatures are, re- 
spectively, RX J1347.5-1145 (Allen, Schmidt & Fabian 
2002 , A2204 (Pointecouteau, Arnaud & Pratt 2005), 
and A1991 (Pratt & Arnaud 2005) respectively. For each 
temperature, we integrate the equations for two values of 
the initial number density: rip = 0.015 and 0.007 cm~^. 
At the initial time the condition ripmpVp / n-xm^V^ = 
(rip/nt)/^^ is assumed to hold. Note that for a NFW 
density profile (Navarro, Frenk & White 1997) the prod- 
uct Uyy^ is independent of position in the central regions. 
We further assume that the gas pressure (oc titT) is equal 
to its initial value. The equations are integrated forward 
for various values of T. 

Figure (1) shows the temperature evolution for for 
velocity independent cross sections, i.e. a = 0. The 
three panels show the temperature evolution for three 
respective values of the initial temperature and initial 
rip of 0.007 (dashed curves) and 0.015 cm~^ (dotted). 
The solid curve in each panel label is obtained with ra- 
diative cooling only. In each panel three different val- 
ues of a have been used. These values, in units of 
10~^^mx/mp cm^, are indicated by the numbers attached 
to the dashed curves. These three values arc also used in 
obtaining the dotted curves (higher values correspond to 



higher curves) . No single value of the cross section is able 
to balance radiative cooling for all temperatures. The 
value of Up clearly affects the time-scale of the evolution, 
with the larger np giving the faster evolution. An inspec- 
tion of this figure allows us to place an upper limit on the 
cross section. The strongest constraint is derived from 
the highest temperature observed clusters of T ^ 15kev 
(Allen, Schmidt & Fabian 2002). According to top panel 
(T^p < 10~^^rnx/mp cm^, consistent with our previous es- 
timate in §4. Exceeding this limit produces temperature 
profiles rising towards the center of all clusters, contrary 
to observations. The limit is insensitive to rip. 

In figure (3), we show the temperature evolution ob- 
tained from the clusters obtained with a^^ = cro(y/Vo)~^ 
where ao = 10~^^(mx/mp) cm^ and Vq = lO^km/s. The 
thin and thick lines correspond to rip = 0.007 and 0.015 
cm~^, respectively. The solid, dotted and dashed are 
for initial temperatures of 16, 8 and 3 kev. The late time 
temperature increase seen for some of the curves is caused 
by the enhancement of heating as a result of virial con- 
traction of the dark matter halo. However, even for initial 
rip = 0.015 cm~^, the temperature increase is minor at 
t = 8 Gyr, corresponding to a redshift of z = 1. Overall, 
we conclude a cross section with a steep dependence on 
energy can balance cooling in clusters of galaxies. 

6. Discussion 

Qin & Wu (2001) have suggested that coUisional en- 
ergy transfer from heavy dark matter particles (rrix/ rrip > 
10^) to protons can balance cooling and thus explain the 
deficit of the cold gas in clusters (Kaastra 2001, Peterson 
2001, Tamura 2001). In the original Qin & Wu scenario 
the cross section is assumed to be a constant independent 
of energy (i.e. the relative velocity between the interact- 
ing proton and DM particle). We have found that a) 
heating of DM with a constant cross section cannot off- 
set cooling for clusters of different virial temperatures, 
and b) any thermal equilibrium obtained with a constant 
(Txp, is unstable. However, an significant upper limit on 
cross sections that are independent of velocity can be ob- 
tained by demanding the temperature not to rise towards 
the center, as seen in observations. Our upper limit on 
a^p versus the mass is represented in the diagram shown 
in figure (4) which also shows previous limits (Wandelt 
et al. 2000). 

Our study incorporates the contraction of dark matter 
halo as it transfers energy to the protons. This contrac- 
tions leads to an increase in the velocity dispersion of 
DM particles and hence to an enhancement of the heat- 
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ing rate. However, for time scales short compared with 
f^jj one can study the thermal stability of the system. 
We have found that for 0-^^ = o-o(V^/Vb)" a stable ther- 
mal equilibrium is possible for a < — 3 (for T > 2 keV 
and metallicity below 0.6 solar). For strong interactions 
the scattering cross-section is expected to have a weak 
dependence on velocity. However, electromagnetic in- 
teractions with a = —4 may fit very well. The DMP 
with a unit charge and mass a few times above 10^ GeV, 
which is just outside the conservative limit on charged 
DMP (De Rujula et al. 1990), seems optimal for sup- 
pressing the cooling flow. Furthermore DMPs of such 
mass if produced in thermal equilibrium would explain 
the present value of fim (Griest & Kamionkowski 1990). 
Another attractive feature of heating by dark matter is 
that it preserves the metal gradient seen in observations 
of some clusters (e.g. Bohringer et al. 2004). This is 
advantageous over mechanical sources of heating such as 
AGN feedback which tend to reduce the metal gradients. 

An alternative possibility for stable thermal balance 
is a DMP with mass comparable to rrip. In this case pro- 
tons and DMPs have comparable energies. Thus no mat- 
ter how large a^^ is, no runaway heating will occur. One 
problem with this scenario is that it leaves a very narrow 
range for rrix- For TOx > 0.4 GeV the allowed values of a^^ 
are too low to achieve thermal balance (Wandelt et al. 
2000). On the other hand DMPs with < 0.2-0.3 GeV 
will typically be less energetic than baryons at the cen- 
tral region. A second problem is that observationally the 
baryon fraction declines towards the center. If the DMPs 
and protons are in thermal balance, then the DMPs with 
mx < l-imp can have steeper density profile only if their 
velocity distribution has significant anisotropy. 

The calculations we have done for protons can be ex- 
tended to other gas particles, such as electrons and he- 
lium ions. Thus for velocity independent interaction we 
find that for helium (TxHc < 1-2 • 10^^''(TOx/mp) cm^, 
while for electrons with Vc = Vp(mp/me)^/^ and = 
1.16np, the limit is a^e < 6 • 10^^^(tox/top) cm^. 

Our analytic treatment is based on the local values 
of physical parameters and therefore can not accurately 
account for the potentially important processes of con- 
vection or the change of central pressure in response to 
the changing temperatures. To describe these effects in 
detail a hydrodynamical simulation is needed. This is 
feasible since the relevant heating function is easy to in- 
corporate in these simulations. 

AN thanks Jerry Ostriker, Joseph Silk and Eti- 
enne Pointecouteau for stimulating discussions. He also 
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Fig. 1. — The temperature evolution of a gas element 
subject to radiative cooling and heating by collisions 
with dark matter particles for constant a^^ {a = 0) and 
Vx/Vp = 0.5. The top, middle and bottom panels cor- 
respond to initial temperatures of 16, 8 and 3 kev, re- 
spectively. The solid curves show are obtained obtained 
with radiative cooling only. Dashed and dotted curves 
are, respectively, for rip = 0.007 and 0.015 cm~^, at the 
initial time. The labels on the dashed curve indicate (Txp 
in units of 10~^^(mx/mp)cm^. 
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Fig. 2. — The upper limit on ctq for Vq = 10^ km/s. The 
solid and the dashed curves correspond to V^/Vp = 1 and 
0.1 respectively. 
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Fig. 3. - The temperature evolution assuming a^^ = 
cra{V/VQ)~'^ The parameter ao is given by ao = 
10^'^^m^/m,p cm~^ and Vb = lO'' km/s. The solid, dot- 
ted and dashed lines correspond to initial temperatures 
of 16, 8 and 3 kev, respectively. Thick and thin lines are 
for initial rip of 0.015 and 0.007 cm~"^, respectively. For 
reference, the look-back time to z=l is 8 Gyrs. 
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